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MONOKINETIC SOLUTIONS TO A SINGULAR VLASOV EQUATION FROM 
A SEMICLASSICAL PERSPECTIVE 

REMI CARLES AND ANNE NOURI 


ABSTRACT. Solutions to a singular one-dimensional Vlasov equation are obtained as the 
semiclassical limit of the Wigner transform associated to a logarithmic Schrodinger equa¬ 
tion. Two frameworks are considered, regarding in particular the initial position density: 
Gaussian initial density, or smooth initial density away from vacuum. For Gaussian initial 
densities, the analysis also yields global solutions to the isothermal Euler system that do 
not enter the frame of regular solutions to hyperbolic systems by R D. Lax. 


1. Introduction and main results 

This paper is concerned by the Cauchy problem for the Vlasov equation 

(1-1) d t f +£d x f - A(d x ln(p))d|/ = 0, f{0,x,£) = f o (x,g), t > 0, (x, £) G R 2 , 

where A ^ 0 and pit,x) = J R f(t,x,£)d£. For A > 0 it arises in plasma physics, e.g. 
for quasineutral plasmas in the core or tokamaks when one focuses on the direction of the 
magnetic lines. There, / denotes the ionic distribution function and the electrons of the 
plasma are assumed adiabatic. 

Due to the derivative of the density p with respect to space in the force term d x ln(p), this 
equation is highly singular. The Cauchy problem can in particular be proven to be well- 
posed for very specific initial data like mono-kinetic distribution functions of the form 

fo(x,0 = Po(x)dx 0 6£ =Vo ( x ) t 

with time-dependent mono-kinetic solutions of the form 

/(f, x, £) = p{t, x)dx <g> 

From a fluid dynamics perspective for A > 0, it is well-known that f(t, x, £) = pit, x)dx ® 
$£=v(t,x) is a distributional solution of (11. Ill if and only if its moments 

pit, x) = J fix, £)d£, and pit, x)v(t, x) = j £/(t, x, f)df 

are solutions of the isothermal Euler system 

j d t p + d x ipv) = 0, 

^ | dtipv) + d x ipv 2 + Xp) = 0. 

This paper addresses the existence of mono-kinetic solutions to ClD, as limits of the 
Wigner transform of solutions to a logarithmic Schrodinger equation introduced in ED 
in the context of wave mechancis. As noticed there, the evolution of initial Gaussian data 
can be computed rather explicitly, a remark which yields our first main result: 
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Theorem 1.1. Let p*,cr o > 0 anduio,po G R. Set 

Po(x) = P*e~ a ° x2 , v 0 (x) = U} 0 x + p 0 , 

and consider the ordinary differential equation 

2A(Jo 


(1.3) 


7 = 


7 


7(0) = 1, 7(0)= wq . 


Suppose that (P) has a solution 7 £ C 2 ([0, T[). Set 

= -^-e~ c ’ o{ - x ~ Pot)2 h{t)2 , v(t,x) = ^jp-x+p 0 . 

7 W 7W 

Then 

P = p{t,x)dx <g> St-yfax) 

is a measure solution to (P > on [0, T[, with p\t=o = podx ® <)j=d 0 . 

In the case A > 0, the solution 7 is globally defined and smooth, so T can be taken 
arbitrarily large. In addition, it safisfies 


7 (t) ~ 2 t\J Acr 0 hit, 7 (t) ~ 2 \/Aero Ini. 

t—foo t—foo 

On the other hand, for A < 0, the solution 7 becomes singular in finite time. 


Remark 1.2. In the case A > 0 the global smooth solutions to the isothermal Euler system 
that are obtained as part of the result, do not enter the frame of regular solutions developed 
by Lax ll24l . Indeed the Lax solutions require bounded initial data, whereas the initial 
velocity in Theorem ll.ll is unbounded. 


Remark 1.3. To our knowledge, the case A < 0 does not correspond to a physical model 
related to (11.11 . However, the logarithmic nonlinear Schrodinger presented in Section [3 
was introduced initially exactly in the case A < 0 (J4]). We will see that the analysis of 
this case requires very little extra effort. 

Our second result deals with more general solutions to d 1.1 b such that p remains bounded 
away from zero. To do so, we do not consider p £ II s (R) the standard Sobolev space, or 
p £ p* + // ,S (R) for some p* > 0, but rather Zhidkov spaces, as introduced in 13 1 11301 1, 
and further analyzed in Ha in the case of Schrodinger equations. Lor s ^ 1, we set 

X S (R) = {/ £ L°°(M), f £ fT s- 1 (R)} . 

Note that being in dimension one, H S (R) C X s (R) holds for all s > 1, and X s (R) is an 
algebra. 

The second result of the paper is the following theorem. 


Theorem 1.4. Let A > 0 and s ^ 2. Suppose that (po^o) € X S (R) x C(K) with 
$0 £ X s (R) and po(x) f Po*for some positive constant po*. 

There are T > 0, p £ C([0, T[; X S (R)), $ £ C([0,T[xR), withd x <5> £ C([0, T[; X S (R)), 
such that 


P = p(t,x)dx ® 5£= dx &(t,x) with (p, $) =(p 0 ,$o), 

t=0 

is a measure solution to o>. 

Moreover, p(t, •) is bounded from below by a positive decreasing function of time. 


Remark 1.5. In (11.Il l, the term d x 
<9? (d x (In / pdC)p) , since f pd£, - 


^ln f pd^jd^p can be considered in a weak sense as 
= p £ C([0, T]; C 1 (R)) and p > 0. 
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Remark 1. 6 . Although the Cauchy problem for the isothermal Euler (11.21) has global in 
time entropy weak solutions (p, it) € L°°(R+ x R ) 2 (see 111 111121 1. we cannot use them 
directly for our purpose. Indeed, the momentum equation in (11.21) is obtained from the 
kinetic equation ( 11 . 11 ) by multiplying ( 11 . Il l by £ and integrating the resulting equation with 
respect to £. This leads to the product 

P d x ln(p). 

Since it is not under a conservative form, it is well known that there is no rigorous way to 
give a sense to this product for a general p £ L°°. It is why we have recourse to regular 
solutions to the isothermal Euler system, thus restricting for general initial bounded data 
to local in time solutions far from vacuum. 

Remark 1.7. For (p,v) with values in ]0, +oo[xR, (11.21 ) is a strictly hyperbolic system. 
Consequently, for (po, t'o) £ (X 2 (R )) 2 with po ^ a for some a > 0, there are T > 0 and 
(p, v) £ (C 1 ([0, T] x R)) solution to the Cauchy problem associated to ( 11.2b and the ini¬ 
tial datum {(>{). vo). This theorem makes an extra connection with a logarithmic Schrodinger 
equation, for which the above system corresponds to a limit system in the semiclassical 
regime. 

Remark 1.8. The assumption on the initial density po is more general than merely po £ 
po* + H s (R) for some s ^ 2, since for instance, po may have different limits as x —> —oo 
and x -A +oo, or, even, no limit at all. 

The plan of the paper is the following. Section[2]recalls the main steps of the derivation 
of the model for A > 0 and related mathematical results. In Section [3] we show how (11.11 ) 
can be obtained formally from a nonlinear Schrodinger equation through the semiclassical 
limit. Section H] establishes Theorem ll.il In Section[5], we prove Theorem l 1 .41 Section[ 6 ] 
adapts the proof of the boundedness from below of the density in the isotropic case by Col 
to the isothermal case. 


2. Derivation of the model and related results 

In this section, we recall the main lines of the derivation of the model with A > 0, used 
for studying fusion plasmas ( fllTl ). The evolution of the ions in the core of such plasmas is 
well described by the Vlasov equation 

dtf + v- V x f + — (-Vx$ + v A B) ■ V„/ = 0, 

where / is the ionic distribution function depending on time, position (in the domain f 1 
of the plasma) and velocity (in R 3 ), Ze and m, are the ion charge and mass respectively. 
The electric potential $ and the magnetic field B should be governed by the Maxwell 
equations. But a finite Larmor radius approximation is derived in the limit of a large 
and uniform external magnetic field. This leads to the following equation for the ionic 
distribution function / in gyro coordinates, 

(2.1) dj + v\\d x J - J° PL {d H *)d v J - (J° PL • V X J = 0. 

Here, the index || (resp. _L) refers to the direction parallel (resp. orthogonal) to the external 
magnetic field. For any vector u = ( m ) £ R 3 , u 1 - denotes the vector (it 2 , —ui, 0). The 
operator J pj is a Bessel operator performing averages on circles of Larmor radius pi, in 
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planes orthogonal to the magnetic field. Since it is not used in this paper, we do not enter 
into more details about it. The electrons move quite more quickly than the ions, so that 
their density n e is given in terms of the electric potential $ by the Maxwell-Boltzmann 
equation 


( 2 . 2 ) 


n e = n 0 e*( M> >. 


Here, e (resp. T e ) is the electronic charge (resp. temperature), and < <t> > is the average 
of the potential on a magnetic field line. Due to the electroneutrality of the plasma, the 
Poisson equation is replaced by the electroneutrality equation 


n e = P, 

where p is the ionic density. The operator J® L induces some regularity in the orthogonal 
direction, but none in the parallel direction. The two-dimensional dynamics in the direction 
perpendicular to the magnetic field is studied in ETI . In order to analyze the difficulty 
coming from the highly singular term J° L (3 X|| ( E , )9„ || /, we restrict to a one-dimensional 
spatial setting, e.g. by considering ionic distribution functions written in the form 

f{t,x,v) = f\\(t,x h v\\)f ± (\vj_\), 

with 

r+oo 

/ f±(\v±\)2ir\v±\d\v±\ = 1. 

Jo 

Then the term f± has no incidence in equation (12.1b and can be factorized. The equation 
that / 1 | should solve is 

d t f\\ +t>|| d x J\\ - A(3 X|| ln( / 0||))d„ || /u =0, t > 0, (at||,V||) e R 2 , 

i.e. the partial differential equation in (11.1b for f\\ (resp. pp, x\\, W||) denoted by / (resp. p , 
x, 0, and A = 


Mathematical results related to (11.1b have been obtained for a system close to equilibrium, 
i.e. in the case where the departure of the electric potential $ from its average along the 
magnetic lines < $ > is small. Equation (12.2b simplifies into 

so that (TO is replaced by 

(2.3) dtf+Cdxf- X(d x p)d^f = 0, t> 0, (x,0gR 2 . 

The Cauchy problem for (12.3b is locally well-posed either for initial analytic data 1221 or 
in Sobolev spaces and satisfying a Penrose stability condition l20l . but is ill-posed in the 
sense of Hadamard for regular initial data in Sobolev spaces and arbitrarily small time 0. 


n e = n 0 1 + 


3. From Schrodinger to Vlasov via Euler 

In this section, we show how a logarithmic Schrodinger equation can be formally related 
to the isothermal Euler system (11.2b . in the semiclassical limit. For e > 0, consider the 
Cauchy problem for the Schrodinger equation 

F i _ 

(3.1) ied t u e + —d^ x u e = X \n(\u E \ 2 )u £ , u e (0 } x) = y / p 0 (x)e l ®°( x ^ s . 
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Following the idea from fUn . any function u e = a £ e l ® e / £ , with (t, x ) a £ (t, x) £ C and 

(t, x ) i—^ , i) el solutions to the quasilinear problem 

(3.2) d t $ £ + (dx * £)2 +Aln(|a E l 2 ) = 0, $ £ (0, a:) = $ 0 (x), 


(3.3) d t a £ + d x ^ e d x a £ + y^xx^ 6 = i^ d La £ , a £ (0,x) = \/po(x) =: a 0 (x), 

is a solution to fl£D - An important remark is that by allowing a e to be complex-valued 
(even though its initial datum is real-valued), one gains a degree of freedom to dispatch 
terms from (13. Il l into (I3.2b - (l3.3b . and the choice introduced by Grenier is much more robust 
than the Madelung transform when semiclassical limit is considered (see 0). 

Determining <t> £ solution to (13.2b turns out to be equivalent to determining if = d x <t> £ 
and a E solution to 

( d t v e + v £ d x v £ + \d x In (|a £ | 2 ) = 0, w £ (0,x) = $' 0 (x), 

(3-4) < n e e „ 

| d t a £ + v £ d x a £ + —d x v £ = i-c% x a £ , a £ (0,x) = a 0 (x). 

Indeed, given (v E , a e ) solution to (13.41) . we can define <f> £ by 


(3.5) $ £ (f,x) = $ 0 (a:) - j Qk(t,x)| 2 + Ain (|a £ (r,x)] 2 )^ dr. 


We check that 

dt (d x $ £ - v £ ) = d x d t $ e - d t v £ = 0, 

so that v e = d x & £ and 4> £ solves (13.21) . and a £ solves (13.31 ). 
Passing formally to the limit e —> 0 in (13.4b . we get the system 


(3.6) 


d t v + vd x v + \d x ln(|a| 2 ) = 0, 
d t a + vd x a + -d x v = 0, O| t=0 


v\t=o = 4>o, 

= t/POi 


which turns out to be the symmetrized version of (11.2b . with p = |a| 2 (see [9] (28J). As 
noticed in the introduction, we then formally obtain a solution to ED by setting 


f(t,x,£) = \a(t, x)| 2 dx ® 

A more direct link from (13.1b to (11.1b is provided by the notion of Wigner measure. The 
Wigner transform of u £ , solution to (14.2b - (14.3b . is defined by (see e.g. ll5l lT6ll26ll ) 

(3.7) W £ (t, x, = f e iyi u £ (t,x- w(t,x+ dy, (t, x, £) e [0, T] x R 2 . 

J M 

Up to the extraction of a subsequence, W £ converges to a non-negative measure on the 
phase space. In general, several limits may exist (see the above references). We will see 
that in the framework of this paper, the limit is unique, and solves ED- 


4. Gaussian initial data 

Let us first study the ordinary differential equation (11.3b . Local existence and unique¬ 
ness of a C 1 solution stem from Cauchy-Lipschitz Theorem. Indeed, the nonlinearity in 
the above equation is locally Lipschitzean away from {7 = 0}. We now address the global 
existence issue. 
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In the case A > 0 , multiplying ( 11 . 3 b by 7 and integrating yields 
( 4 . 1 ) (i) 2 = w o + 4A<7o In 7. 

This readily shows that 7 is bounded from below away from zero, so the flow is global. 

In the case A < 0 , suppose that 7 is bounded away from zero, i.e. there is S > 0 such that 
7(f) ^ S. Then ( 11.3b yields 

.. 2Auo 

hence 7 (t) ^ pp t 2 + loq t + 1, and a contradiction for t sufficiently large. Now suppose 
that 7 £ C 2 ( 0 ,00) with 7 > 0 : from the above argument, there exists a sequence t n along 
which 7 (t n ) —> 0 + . From (11.3b . 7(i n ) —> —00, hence a contradiction. 

To prove Theorem l 1.11 we start from a semi-classically scaled logarithmic Schrodinger 
equation, 


(4.2) 

together with the initial value 


ied t u e + —d 2 xx u e = A ln(|M E | 2 )-u E 


o -a 0 x 2 /2 iu> 0 x 2 /(2e) ipox/e 


(4.3) u E (0,x) = 

Such an initial datum does not fit into the framework of |j2]|, since it goes to zero at infinity. 
However, as noted in j4j, for fixed e > 0, the solution to (14.2b with such an initial datum 
can be computed rather explicitly. Indeed, bearing in mind the propagation of coherent 
states in the semi-classical regime ( |[T9l , see also |fT3l ) it is consistent to look for a solution 
of the form 

( 44 ) u e {t,x ) = 

with b £ , fl E £ C, and q,p, S £ R. Plugging this ansatz into (13.1b . we find: 

ieti - iefr <KX ^ q) \ £ + ieQ £ q{x - q)b e - p(x - q)b e + pqb e - y b e - Sb e - e 2 y b e 


■ £ 2 (H 


c\2 


(x-q) 2 


b E - ipen e (x — q)b E = A b E In |6 E | 2 - A b E {x — q) 2 Re fl E . 


Cancelling the polynomial in x — q, sufficient conditions for the previous equation to hold 
are 


(4.5) 


q= P , p = 0, 5 = — y = y, 


ieb E = y f Yb e + A b E In \b‘ 


e\2 


(4.6) 

(4.7) ietl E = £ 2 (H e ) 2 + 2AReH E . 

Equation (14.5b corresponds to the classical Hamiltonian flow in the absence of external 
force, and the associated classical action. We compute exactly 

P(t)=p 0 , q(t)=p 0 t , S'(f) = yf. 

Since we are eventually interested only in the modulus of b E , we infer from (14.6b : 


hence 

(4.8) 


y| 6 E | 2 = 2Reb E b E = s\b E \ 2 ImH e , 


\b E (t)\ 2 = p^ e e /o Imfie (^7 
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For fixed e > 0 and Re fl E (0) > 0, (14.71 ) has a unique, global solution fl E £ C(R), whose 
large time behavior depends on the sign of A. Indeed, in a similar fashion as in li25j . we 
seek Q. £ of the form 


n £ = - 


Then ( 14.7b becomes 


i - 

er] E 


-e 2A 7 f 

rf — — ?? Im —. 
e rj e 

Introducing the scaled polar decomposition r] s = 7 E e lEI 
(4.9) 


is given by 


n £ = a E - 1 -^. 

£ 7 e 


and the above equation reads: 

7 e — e 2 7 e (d e ) 2 = 2X'y s 
The second equation yields 


ci E 7 E + 2a EA f e = 0. 


dt 


(( 7 e ) 2 d e )=0. 


In view of (14.91 ), we have a degree of freedom to set 7 £ (0). Setting 7 £ (0) = 1, we have 

a e (0) = Refl E (0) = opi 7 £ (0) = ^eImfl E (0) = to o- 
Therefore, we have 

(4.10) ( 7 E ) 2 d E = cr 0 . 

The equation for boils down to 

" 2 -'o 


7 = 


£ 2 <t 2 2Acto 


( 7 e)3 




along with the initial data that we recall: 

7 e (0) = l, 7 e (0 )= Wo . 

Local existence and uniqueness of a C 1 solution stem from Cauchy-Lipschitz Theorem. 
Indeed, the nonlinearity in the above equation is locally Lipschitzean away from { 7 E = 0}. 
We now address the global existence issue. Multiplying the above equation by 7 E and 
integrating in time, we infer: 


(7 e ) =4A(7oln7 E — 


£ a. 


(7 E ) 


0 , 22 , 2 
+ £ °n + w n- 


This shows that for fixed e > 0, 7 e remains bounded away from zero, for if we had 
7 e (t n ) —> 0 for some sequence t n , then the above right hand side would become negative 
for n large enough, hence a contradiction. Thus, for fixed e > 0, (14.71) has a unique, global 
solution fl E £ C 1 (R). 

In view of (14.101 ). (14.91 ) also reads 


(4.11) 


n e = 


do 

(7 e ) 2 


i 7 
£ 7 e 


Given (14.81 ) and (14. 1 lb . the Wigner transform (13.7b of u £ is given by 


n _ ^n^-w 1 ) 2 

(4.12) W e {t,x,0 = —r^e h'W) 2 

7 £ (t) 


e iv(.i~ (x-Pot)-po) d 
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On every time interval such that 7 is bounded away from zero, the Gronwall lemma shows 
that 


7 e - 7 = 0(e), 7 e - 7 = e>(e). 

Consequently, when s —> 0, the Wigner transform W e of u e weakly converges to the 
bounded measure 

p* _fo( x — po*) 2 


P _ ■t'U '-J 

K*,dx,dO = ^e w*» dx®6^ x _ pot)+po . 


(4.13) 

Straightforward computations show that 

/ P * gp^-po *) 2 / y(t) \ 

(p(t,x),v(t,x)) := (^y e , ^)( X ~ Pot ' ) +Po ) 

is a solution to the isothermal Euler system (11.2b on the time interval of the existence of 7. 

To conclude, we consider the large time behavior of 7 for A > 0. If 7 was bounded 
from above, 7 ^ M, then ( 11.31 ) would yield 

•• ^ 2A<t 0 

7 ? Tir >0 ' 

hence a contradiction. Therefore, for t sufficiently large, 7 (f) ^ 1 and 7 (f) > 0. Since 
7 ^ 0 , integration then shows 

7 (t) > i(t)(t ~t) + 1 —*• 00, 

t—> OO 

and 7 > 0 for t sufficiently large. The asymptotic behavior announced in Theorem 11.41 
then follows by integrating the identity 

d'y 

V^o + 4A(Jo In 7 

to obtain the asymptotic behavior of 7 . 

d'y 


= dt. 


t = 


y/ojQ + 4A<7o In 7 2 A<to 
where we have changed the variable as 

V = \j 4Afj 0 ln 7 + Wq. 

Recall that the Dawson function, defined by 


— [ e (y2 "“o )/(4A<To) dy, 

Wo J 


satisfies 


F(x) = e 
F(x) 


x —>+00 


e y dy 

1 

2x’ 


(see e.g. ID), we infer that 


Since 


we conclude 


7 (t) ~ 2 t\J Ado Inf- 

t—> + 00 

7 = \]ul + 4Acr 0 In 7 , 


7(f) ~ 2y/Aero Inf. 
£—>■+00 
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In the case A < 0, suppose that 7 is bounded away from zero, 7 (f) ^ 5 > 0. Then (11.3b 
yields 

.. 2A(Jo 

1 ^ S’ 

hence 7 (t) ^ t 2 + Uq t + 1, and a contradiction for t sufficiently large. Now suppose 
that 7 G C 2 (0, 00 ) with 7 > 0: from the above argument, there exists a sequence t n along 
which 7 (t n ) —>■ 0 + . From (11.3b . ■y(tn) —> — 00 , hence a contradiction. 

This ends the proof of Theorem ll.il 


5. WKB ANALYSIS 

In this section, we justify the formal approach presented in Section[3]in the framework 
of Theorem ll.4l thus proving this result. 


5.1. Constructing the solution. Note that for fixed e > 0, the Cauchy problem for (13.1b 
has been considered in J 8 ] (see also 13 Section 9.1]), for initial data in the class 

W = {/ G H\R), \f(x)\ 2 |In \f{x)\ 2 \dx < 00 1. 

This class is not compatible with the assumption p(x) ^ po* > 0 from Theorem 11.11 
which is equivalent to |u e ( 0 , x)\ 2 = |o,q( at) | 2 ^ po* > 0 in the approach that we follow. 
Therefore, we choose to rather work in Zhidkov spaces X s (K). The system ( 13.4b has a 
unique smooth solution as stated in the following proposition, which includes the case 

£ = 0 . 


Proposition 5.1. Let s > 5/2 and A > 0. Suppose that po, $q G X s (R), with 

Po{x) > Po* > 0. 

Then there exists T independent of s > 5/2 and £ € [0,1], and a unique solution 
(, a E ,v e ) G C([0, T]; X s x X s ) to (lT4b. 


Proof. This result is a rather direct consequence of [(21 Proposition 2.1], whose proof we 
recall the main idea. Separate real and imaginary parts of a e , a e = a\+ ia |, and introduce 


u 


£ 



u 0 = 0 

V^O / 


L = I dl x 


—d 2 

T.T 

0 

0 


and Al(u) 


We now have the system: 


v 0 

0 v 

2 A a 1 2 A a, 2 

a l + a 2 Ct l - l -a 2 



(5.1) 


d t u £ + A(u e )^u £ = ^Lu £ ; u [ t=0 = u 0 . 


Since po is bounded away from zero, its square root is also in X s , so that Uo G X S (K) 3 . 
The matrix A is symmetrized by the matrix 


S = 


h 

0 


0 

a^+al 

4A 
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which is symmetric positive if and only if a\ + a| > 0, that is, so long as no vacuum 
appears. By assumption, 


(«D 2 + (a e 2 ) 2 


^ Po* ^ 0- 

t=o 


Then the main idea is that the operator L is skew-symmetric, and so does not appear in 
/Abased energy estimates. Standard tame estimates (see e.g. 12711291 ) do not involve the 
L 2 norm of u e , and so the only aspect remaining is that L°°-estimates can be obtained 
rather directly. So long as, say, 


(5.2) (a\(t, x)) 2 + (a|(f, x)) 2 ^ Vx G K, 

we have: 


||u e (/)||l=o < ||u 0 ||l=o + f \\A(u 6 (T))d x u s (T)\\ L ~dT + f \\d 2 xx a £ {T)\\ L ^dT 

J 0 JO 

< ||u 0 ||l~ +c f ||u £ (r)|| L =o||a a; u e (r)|| L oo C /T+ f \\d 2 x a e (T)\\ L °*dT 

J 0 JO 


sdr+ [ ||u e (T)\\x»dT, 

Jo 

where we have used Sobolev embedding under the assumption s > 5/2. Indeed, by def¬ 
inition, X s C L°°, and for u e G X s , d x u e G 7T S_1 C L°°, provided that s > 3/2, and 
similarly, d xx \i £ G L°° for s > 5/2. 

Now we set P = (I — d xx ) < ' s ~ 1 ^ 2 d x , so that ||/||x s ~ ||/||l°° + ||P/||i 2 , and denote by 


^ H u oIIl“ + c 




?x 



f{x)g{x)dx, 


the scalar product in L 2 . Since L is skew-symmetric and S is real-valued, 

p(SPu E (i),Pu E (i)> 

= {(dtS)Pu E {t),Pu E {t)) + 2Re(Sd t Pu £ (t),Pu £ (t)) 

= {(d t S)Pu £ (t),Pu £ {t)) + eRe(SLPu £ (t) : Pu £ {t)) 

— 2 Re ^ SP^A(u e (t))d x u e (t , Pu E (/)^ . 

So long as (15.21 ) holds, we have the following set of estimates. First, 

((<9 t S)Pu E (f),Pu E (/)) < \\d t S\\ L ~>\\Pu £ (t)\\ 2 L2 

<C(||u E (/)||^)||a t u E (/)||^||u E (/)||^. 

Directly from (15.1b . we have: 

II^u e (/)||l=° < C (||u E (/)||ioo) \\d x n £ m^ + \\d 2 x a £ (t)\\ L ~ 
^C(\\u £ (t)\\ x ,)\\n £ (t)\\ X s. 


Since SL is skew-symmetric, we have 

Re (SLPu £ (t), Pu e (t)) = 0, 


which prevents any loss of regularity in the estimates. For the quasi-linear term involving 
the matrix A, we note that since SA is symmetric, commutator estimates (see ||23l ) yield: 

(5P(A(u £ )a iU E ),Pu E (t)) < C , (||u E (t)|| L ») ||Pu e (f)|||2||5 x u E (i)||i,=o 

<c(||u E (/)|U 8 )l|Pu e WII 2 L2- 
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Finally, we have: 

^(SPn%t),Pn%t))^C(\\n%t)\\ x ,)\\n%t)r x ,. 

This estimate, along with the L°°-estimate, shows that on a sufficiently small time interval 
[0, T\, with T > 0 independent of e £ [0,1], (15.2b holds, hence the existence of a unique 
solution. The fact that the local existence time does not depend on s > 5/2 follows from 
the continuation principle based on Moser’s calculus and tame estimates (see e.g. EH 
Section 2.2] or ll29l Section 16.1]). □ 


Corollary 5.2. Under the assumptions of Proposition 1377] if we suppose in addition that 
$q £ L 2 (R), then ( 13.1b has a unique solution u e £ C([0, T]; X S (R)), where T is given by 
Proposition 15.71 


Proof. From Pronosition l5.il (13.41) has a solution ( v e , a e ) £ C([0, T]; X s x X s ). Plugging 
this information into (13.4b . we infer 


< KIL’+ f \\v E (T)\\ L ~\\d x v s (T)\\ L2 dT 

Jo 

\\d x a E (T)\\ L 2 dT. 


+ C 


r 

1 

L 

a e (r) 


Therefore, v e £ C([0, T]; L 2 ), and 4> £ , stemming from v e via the formula (13.5b . satisfies 
<f> e £ C([0, T]\ X s+1 ). The existence part follows readily, since X S (R) is an algebra, by 
setting u e = a E e l * ! e . 

For the uniqueness property, consider two such solutions u e , vr £ C([0, T]; X s ), and 
set w e = u e — u E . It satisfies 


£ 2 

(5.3) ied t w E + — d 2 x w E = A (ln(|u £ | 2 )u £ — ln(|u e | 2 )u e ), w E t=0 = 0. 

Recall the pointwise estimate from [|8) (see also El Lemma 9.3.5]), 

|lm (ln(|u e | 2 K - ln(|u e | 2 )u e ) ( u e - u E )\ < 4| u e - u e \ 2 . 

Multiply (15.3b by w e , integrate on an interval / = [M_, M + \, and take the imaginary part. 
This yields, along with the above estimate, 

/ \ wS (t,x)\ 2 dx + E — !m J nfd 2 x w e ^ AX Jjw E (t,x)\ 2 dx. 

We have, by integration by parts, 


Im 


w e d 2 w e 


Im w e (t, M+)d x w E (f, M+) — Imw e (f, M-)d x w e (t, M_). 


Since w E £ C([0, T]; X 1 ), we can choose sequences M” —> ±oo along which the above 
term goes to zero, and the Gronwall lemma implies ||(i )||^2 =0. □ 


5.2. Asymptotic expansion. Proposition [5T| with e = 0 yields the existence of a unique 
solution (v, a) £ C([0, T]; (X S (R)) 2 ) to (13.6b As a direct consequence of Proposition l5.il 
and El Proposition 3.1], we have: 

Proposition 5.3. Under the assumptions of Proposition 13771 there exists C independent of 
e £ [0,1] such that 

\\dx — $) IIl“([o,t] ; x s - 2 ) + ll aE — a llL“([o,T] ; x: a - 2 ) ^ Ce. 
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5.3. Convergence of the Wigner transform. For u e = with (<1> £ . a e ) solution 

to (13.2b — (13.3b . or equivalently (d x <& e , a e ) solution to (13.41) . the Wigner transform defined 
in (13.71) is equal to 


(5.4) W £ (t, x, £) = J e^ v a £ (t, x - | y)W (t, x + |y) e^^^dy, 

where 

<p E (t, x, y ) = (t, x - |y) - (t, x + |y) . 

Theorem 5.4. When e —> 0, the Wigner transform W e of u e weakly converges to the 
bounded measure 

y(t,dx,d£) = |a(i,x)| 2 dx® S^ dx$(ttX) , 
where (d x <f>, a) is a solution of (13.61) . Moreover, y is a solution to (P) . 

Proof. In view of Proposition [53] a s = a + r £ and d x <f> s = d x 'P + rf with 


ll r all-£/’ o ([0,T];X' i - 2 ) + lku||z,°°([0,T];X'— 2 ) ^ Ce. 


Therefore, 

W £ (t,x,0 = j (t, x- £ -y)d (t, x + |y) dy + R\ + R s 2 + i?|, 

where 

Rj(t,x ,£) = J e lv{i ~ d ^ ( ' t ’ x)) r £ (t,x,y)dy, 1 < j < 3 , 

rf (t,x, y) = a £ (i, x - |y) F (i, x + |y) W,*))/e - l) , 

r e 2 (t,x,y) = d (t,x+^yj r e a (t,x- |y) + a (f,x - |y) rf (f,x + , 

^I(*i ^ 2/) = r a (*, x-^yjrf (t, x + |y) . 

Propositions [53] and [53] yield, along with Taylor’s formula for the term r \, 

II^||l~([o,t]xr 2 ) < Ce, 1 < j ^ 3. 


Consequently, W e tends to |a| 2 <ix 0 <5j=g x $ in A4f,([0, T] x R 2 ) when e tends to zero. 
Moreover, denote by (•, •) the duality between bounded measures on [0, T] x R 2 and contin¬ 
uous functions with compact support in [0, T] x R 2 . For any test function a(t, x, £) £ C 1 
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with compact support in [0, T] x R 2 , it holds 


(/a d t a + £d x a - A (d x ln|a| 2 ) d$a) 

, x)\ 2 (■ d t a(t,x,d x $(t,x )) + d x $(t,x)d x a(t,x,d x $(t,x))) dxdt 



— Ay |a(i,x )| 2 ((Sc ln|a| 2 ) d^a(t,x,d x ^(t, x))) dxdt 
/ |a | 2 (dt (ac(t, x, v(t, x))) — drvd^a/t, x, v(t , x))) dxdt 


+ / |a | 2 (d x (v a(t, x , i>(i, x))) — a(t, x, v) — vd x v d^a{t, x, v)) dxdt 


— Ay |a| 2 ((9 X In|a,| 2 ) d^a(t, x, v(t, x))) dxdt 

— j a(t, x, v) (d t \a\ 2 + vd x \a\ 2 + \a\ 2 d x v) dxdt 

— / \a\ 2 d^a(t, x, v ) (d t v + vd x v + \d x ln|a| 2 ) dxdt. 


This is zero, in view of (13.6b . since 


d t \a\ 2 + d x $>d x \a\ 2 + \a\ 2 d xx $ = 2Rea (d t a + d x <bd x a + |9 2 X T>^ = 0. 


Therefore, any solution to (13.6b yields a solution to dl.lt . 


□ 


6. The bound from below of the density 


It remains to prove that the density is bounded from below to complete the proof of 
Theorem [T4l The result from IfTOll yields: 

Proposition 6.1. Under the assumptions of Theorem \1.4\ the density p solution to dl.2t 
satisfies 


pif) > i P + Ct ' *G[ 0 ,+oo[, 


for some constant C > 0. 

Proof In flOl , the proof is presented for the pressure law p(p) = p in ( 1 1 .2b replaced by the 
isentropic one, p(p) = (P , with 7 > 1. We simply perform the standard modification to 
adapt the proof to the isothermal case (see e.g. imiT4l ). Note that Proposition 16. II is also 
a consequence of Col Remark 3.3]: we sketch the proof in the isothermal case for A = 1 
for the convenience of the reader. 

First, the isothermal Euler equation in Lagrangian coordinates reads (see e.g. Ifl4l ) 


( 6 . 1 ) 



T 


with t = 1/p. Setting 


s = u — In r, r = u + In r, 


we find 


<9+s = d-r = 0 , 



T 


where 
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Following flOl. introduce 

a = s x , /3 = r x . 

By assumption, we know that 
exists M > 0 such that 

u, t, a and (3 are uniformly bounded at time t 


sup(a(0,a;),/3(0,x)) < M. 


x€M 

We check 

d+a = ^a(/3 - a), 

( 6 . 2 ) 

(6.3) 

<9-/3 = ^r/3(a - 0). 

As in m3. 


(6.4) 

sup (a(£, x), /3(£, x)) < M, 


(t,x)G[0,T)xR 


where T > 0 is such that (16.1b has a C 1 solution on [0, T). (16.4b is established by contradic¬ 
tion: assume for instance that a(£o, ^o) = M. Because the wave speed - is bounded from 
above, we can consider the characteristic triangle with vertex (f (l , xq) and lower boundary 
at time t = 0, denoted by Cl. Let fi(^ to) be the first time such that a(ti, Xi) = M or 
0{ti,Xi) = M in fi. Assume a(fi, a;i) = M for instance (the other case is similar), and 
let Cl\ denote the characteristic triangle with vertex (t-\. X\): 

sup (a(t,x), /3(t,x)) < M. 

(t,x)€£ 2i, t<ti 

Since a(ti,x±) = M, there exists t 2 £ [0,£i) such that 

inf a(t,x) > 0 . 

(£,:r)6f2i, t^t2 

Let t 2 < t < t\. On [t 2 , t ], (16.2b yields 

d+a ^ K(Ma — a 2 ). 


Integrating along some charateristic between £2 and t, we find 


— In a{t) 

M M — a(t) 


< J_ l n a ( t2 

^ M M- a(t 2 ) 


+ K (t — £ 2 ) • 


Letting t —> t\ then leads to a contradiction. 

We infer that u x is bounded from above, hence r t is bounded from above, in view of 
dO : t grows at most linearly in time, and since p — 1 /r, the proposition follows. 

□ 
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